X-ray observation of micro-failures in granular piles approaching an avalanche 
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An X-ray imaging technique is used to probe the stabihty of 3-dimensional granular packs in a 
slowly rotating drum. Well before the surface reaches the avalanche angle, we observe intermittent 
plastic events associated with collective rearrangements of the grains located in the vicinity of the free 
surface. The energy released by these discrete events grows as the system approaches the avalanche 
threshold. By testing various preparation methods, we show that the pre-avalanche dynamics is not 
solely controlled by the difference between the free surface inclination and the avalanche angle. As a 
consequence, the measure of the pre-avalanche dynamics is unlikely to serve as a tool for predicting 
macroscopic avalanches. 

PACS numbers: 45.70.Cc, 45.70.Ht, 81.05.Rm 



When a granular pack is submitted to a slowly vary- 
ing stress its apparent response consists in intermittent 
bursts of plasticity during which large irreversible defor- 
mations take place 0, |3| • Between these discrete events 
the system behaves as a rigid body. Formation of shear- 
bands in granular systems in triaxial tests the jam- 
ming and unjamming of an hourglass 0] or sequences of 
avalanches down a sandy slope are common examples of 
this intermittent behavior. Because of its relevance to 
geophysics and the possibility of direct observations, the 
case of avalanches have been studied in great detail since 
the early work of Bagnold 0. It was found that the 
intermittent regime occurs between two limiting angles 
of the heap surface: the repose and the avalanche angle 
which depend on the geometry as well as the characteris- 
tics of the material (shape and surface properties of the 
grains) ^ . Although the main features of the avalanche 
flow are rather well understood 8, 9] , one important 
question remains to date unsolved: what is the nature 
of the microscopic process by which the flow is triggered 
and later stopped during a single avalanche? 

To address this question several experimental and nu- 
merical studies have recently focused on micro-plasticity 
[lO, 11]. Indeed it appears that even in the absence of vis- 
ible flow, any modification of the external force applied 
to a granular pack induces some micro-displacements 
which allow the system to mechanically adapt to the 
new constraints. By analogy with earthquakes, these 
events could be seen as precursors of the macroscopic 
failure of the pack. Although they dissipate only a small 
fraction of the energy as compared to the macroscopic 
avalanches, they could impact the stress distribution 
within the static pack therefore controlling its overall 
stability. 

In the present work, we use an X-ray imaging tech- 
nique to probe the energy and spatial characteristics of 
these precursors in 3D granular systems slowly tilted to- 



wards the avalanche angle. Experiments are done in a 
cylindrical acrylic drum 55 mm in diameter with a depth 
of 30 mm. The drum is half filled with various types of 
grains. Table and mounted with its rotation axis par- 
allel to the X-ray beam. To prevent a bulk rotation of 
the sample, a disordered monolayer of grains is glued to 
the inner surface of the cylinder. The drum is co-axially 
mounted inside the central aperture of a precision rota- 
tion stage which allows for very smooth rotation and a 
maximal rotation speed of 0.5°s^^. A micro-focus X-ray 
source is placed 230 cm from the drum and provides a 
weakly divergent beam which crosses the stage aperture, 
producing a radiograph on a scintillator-coupled CCD 
camera 12], figure^a). For a monochromatic beam, the 
transmitted intensity of the incident beam at a position 
\x, z] reads: 



I[x, z] — A(t) ■ Iq[x, z] ■ a.c[x, z] ■ exp 



Lg[a:,2:] 



(1) 



In this expression, A(t) ■ 1q[x,z] is the incident intensity 
map without sample {A{t) being the temporal fluctua- 
tions of the X-ray source intensity), ac[a;,2:] corresponds 
to the absorption of the empty container, and Xg is the 
absorption length of the grain material. The total length 
of the grain material which projects through each posi- 
tion [x,z] on the drum is denoted, Lg[a;, z]. 

To prepare the grain pack in a reproducible way, the 
drum is initially rotated 360° clockwise, leaving the pack 
surface at the avalanche angle —9a- The drum is then ro- 
tated incrementally counter-clockwise. To probe whether 
a displacement has occurred over any given increment SO 
(0.5°), we use the following protocol: at a radiograph of 
the pack is taken, Ii [x, z]{6). The drum is further rotated 
to the angle 9 + 59, stopped for 5 s, then rotated back 
to 9 whereupon another radiograph l2[a::, z](0) is taken. 
Finally, the drum rotates back to 9 + 69 in preparation 
for the next increment. The change in lig[x, z], by direct 
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FIG. 1: Measurement of the dynamics induced by a tilt in- 
crement for the system I (see Table HJ: (a) radiography of 
the drum, (b) picture presenting the variations of the crossed 
grains length ALg, at an angle 20°, (c) measurements of the 
quadratic fluctuations as a function of the drum angular po- 
sition 0. The doted line corresponds to a single run, whereas 
the inset (semi-log scale) shows the average over 15 runs. The 
data have been rescaled by the average of Q in the avalanche 
regime. 





Material 


Shape 


Size 


ea 


I 


Pasta 


cylinder 


2.2mm 


30 


II 


Carbon 


sphere 


1.1mm 


27 


III 


Carbon 


sphere 


0.8mm 


26 


IV 


Carbon 


sphere 


0.4mm 


22 



TABLE I: Different types of grains used in the experiments. 
9a is defined as the average inclination angle of the sur- 
face in the avalanche regime. The size corresponds to the 
diameter for the spheres. The cylinder is 2.2 ± 0.2mm in 
height and diameter. Glassy carbon beads were obtained from 
Hochtemperatur-Werkstoffe Gmbh, pasta was De Cecco Acini 
di pepe. 



comparison between Ii[a;, z](0) and l2[a;,2;](0) quantita- 
tively characterizes the pack dynamics: 



AL,(0) - A, 



In 



h{0) 



In 



hi0) 



(2) 



The average in eq. El corrects for temporal fluctua- 
tions in the X-ray source. A grey-scale representation of 
ALg[a::, z] is presented in figure ^b). This image demon- 
strates the existence of movements in the vicinity of the 



free surface. A closer examination reveals that these mo- 
tions occur in the form of clusters of a few grains, moving 
over a fraction of the grain size. 

Several measurements may be extracted from 
ALg[x,z]. The quadratic average of ALg{9), denoted 
Q{0), is a probe of the amplitude of small relative grain 
movements (a fraction of the grain size). Since this 
average also integrates the detector noise, we evaluate 
the contribution of the latter, denoted Qo, by measuring 
the quadratic noise on successive pictures without 
rotation and taken under the same conditions as in the 
experiment. We subtract this value so that Q{0) only 
reflects grains motions: 



Q{0) = {ALg{0)') ~ Qo 



(3) 



We can also extract, for the same measurements, the 
displacement (Axb, Azh) of the barycenter of the pack, 
for every angular increment, 60: 



( Axb{0) 
[ Azk{0) 



= K 



E 

{x,z} 



A-Lg{e)[x,z] (4) 



The sum is performed over all the positions [x, z]. K is 
chosen in order to express the distances Axh and Azf, 
in drum radius unit. Though expressed as a length, Az 
may also be interpreted as an energy variation. 

Figure^c) displays the evolution of the quadratic fluc- 
tuations Q{0)- Its spikes are a direct signature of the 
existence of intermittency in the micro-movements of the 
grains, well before the surface reaches the avalanche an- 
gle. This observation is consistent with the 2D numerical 
simulations of Staron et al. ^3 ■ The discrete events be- 
come detectable when the surface inclination approaches 
the horizontal. We thereafter observe a progressive in- 
crease of the magnitude of Q{0)^ but no significant evo- 
lution of the intermittency frequency. 

The displacements of the pack barycenter, on differ- 
ent systems, and with different grains sizes, are dis- 
played in figure El Three regimes may be distinguished. 
For — < < 0, we measure small values of Az^, 
whereas Ax^ remains insignificant. This indicates that 
a slight compaction occurs in the sample as expected 
when a granular system is mechanically perturbed. For 
< 6' < 0a, we observe a strong increase on both Axb 
and Azh- This evolution is associated with a progres- 
sive increase of the mean displacement direction, given 
by tan~-'^(Aa:h/Azf,), that eventually aligns with the free 
surface, figure (2 inset. As the inclination reaches the 
avalanche angle 0a^ we observe a sudden increase of the 
dynamics; beyond 0a, the flow occurs along the surface of 
the pack ((tan~-^(Aa:b/Azb)) = Oa). The values of Axb 
and Azft simply correspond to the average grain flow nec- 
essary to compensate for the rotation of the drum. 
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FIG. 2: Normalized displacement of the pack bary center for 
four different systems. The length unit is set as the drum 
radius. (I) is the average over 15 runs, whereas (II), (III) 
and (IV) are averaged over 25 runs. The inset describes the 
angular direction of the displacement. 



Surprisingly, the measurement in the pre-avalanche 
regime seems to be independent of the grain type, and 
in particular of the grain size. This observation can be 
used to discriminate between two possible interpretations 
(see figure |3l : (a) If the dynamics were due to a purely 
superficial effect, such as grains rolling, the thickness of 
the moving layer and the magnitude of the displacement 
would be controlled by the grain size. Hence, the aver- 
age displacement in the pre-avalanche regime, expressed 
in units of drum radius, should scale as the square of the 
grain size. This is inconsistent with our observations, 
(b) In contrast, the independence of the dynamics with 
grain size strongly suggests that the surface flow pene- 
trates onto a depth d controlled by the drum geometry 
(diameter and/or thickness). The pre-avalanche regime 
thus corresponds to a macroscopic instability of the pack 
structure, though preferentially occurring in the vicinity 
of the surface. 

In considering figure |21 it is tempting to describe this 
instability as being controlled by the surface angle. To 
probe this hypothesis, we explore the effect of other pa- 
rameters - such as the preparation - on the pre-avalanche 
dynamics. The following sequence of tilt is thus applied 




superficial displacement deep reorganization 

FIG. 3: Illustration of the superficial and deep displacements. 
Case 1 (left): the displacement involves a superficial layer 
which thickness scales with the grain size r, moving over a 
distance of the order of r. The barycentric displacement Ax 
then scales like r^. 

Case 2 (right): the failure corresponds to a minute deforma- 
tion of a layer whose thickness d is controlled by the drum size: 
in this case, the barycentric displacement is independent of 
the grain size. 



to the pack before starting the measurement: the pack 
is first rotated 360° clockwise, then slowly tilted back- 
wards until the surface inclination reaches a prescribed 
angle ^o, and finally returned to a fixed value —10°, fig- 
ure ^a). From that point, the dynamics are recorded 
as the drum is rotated counter-clockwise using the usual 
protocol. 

The figure Efa) shows the evolution of the quadratic 
fluctuations Q{9) (equation O obtained for three differ- 
ent values of Oq, which are compared with the reference 
curve obtained with data from figure ^inset. Interest- 
ingly, the dynamics are, in all cases, negligible until the 
drum reaches the angle ^o- Beyond this value, the activ- 
ity increases rapidly and shows similar features as in the 
non pre-loaded sample. This result demonstrates that 
the pre-avalanche dynamics does not depend solely on 
the surface inclination. The absence of detectable plastic 
events from —10 to 9q shows that the successive static 
configurations that the system experiences during a tilt 
experiment are stable over a large range of surface incli- 
nations. 

To probe in more detail the range of stability of 
each structure we perform a complementary experiment 
depicted in Figure El^b). As previously shown, the 
system is rotated 360° clockwise, then tilted backwards 
to a prescribed angle Oq < 0a- The dynamics is then 
recorded as the pack is rotated clockwise towards the 
avalanche angle. We observe that the threshold angle 
for which detectable motion occurs, and consequently 
the range of mechanical stability of the pack, depend on 
the preparation angle 6*0. For 6*0 of the order (or less) 
than 0, the pack remains stable until we reach —9a- 
For larger values of ^o, the onset of internal dynamics 
occurs after a rotation of roughly 30°, regardless of 
9o. These results indicate that each structural state 
of the pack is associated with a range of inclination 
angle within which it is mechanically stable (within the 
resolution of our measurements). The width of this 
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In contrast, the avalanche is a dynamical instability that 
primarily involves the dissipative properties of the mate- 
rial. An avalanche is triggered when the pile cannot dis- 
sipate the inertial energy produced by a any small failure 



sipate tne 



We have presented here a new experimental technique 
suitable for the study of 3D confined granular dynamics, 
that allowed us to measure minute displacements of the 
structure. In contrast to other techniques (DWS ca- 
pacitive measurements ^^), our method allows for local 
and quantitative measurements of the deformation. We 
are currently using X-ray computerized tomography in 
order to obtain 3D information on these plastic events 
[l7| |. Beyond the displacement field, this method might 
enable us to directly probe the statistical properties of 
the contact network, which might control the plastic re- 
sponse of the material. 

We wish to thank Tomaso Aste for stimulating discus- 
sions, and Tim Sawkins for his helpful technical contri- 
butions. 
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FIG. 4: Role of the preparation (grain type I) (a) Effect of a 
pre-load until an angle 9o on the dynamics, (b) Evaluation of 
the range of stability for different preparations. All the curves 
correspond to the average of four similar runs (except 15 for 
the reference curve) . For a better visibility, graphs have been 
shifted. 

stability range barely depends on the pre-charge an- 
gle, but the latter controls the value of its limiting angles. 

From the observed pre-avalanche dynamics, one might 
be tempted to extrapolate information about the trigger- 
ing of macroscopic avalanche flows such as the avalanche 
angle. However, we find that grains with differ- 
ent avalanche angles exhibit quantitatively similar pre- 
avalanche behaviors, figure |21 Moreover, we show that 
the preparation and history of the pile have a major im- 
pact on the statistics of the micro-events, but no signifi- 
cant effect on the avalanche angle As a consequence, 
it seems unlikely that the measure of the pre-avalanche 
dynamics may serve as a tool for a robust prediction of 
the macroscopic avalanche. 

These observations might reflect the fact that the two 
instabilities - pre-avalanche failures and avalanche trig- 
gering - are qualitatively different. The discrete micro- 
failures we have observed correspond to the mechanical 
destabilization of a static structure induced by an incre- 
mental tilt. The energy released by each event is instan- 
taneously dissipated. The instability is thus only con- 
trolled by the static properties of the granular structure. 
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